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ABSTRACT 

The Cosmic InfraRed Background (CIRB) is the sum total of the redshifted and 
reprocessed short wavelength radiation from the era of galaxy formation, and hence 
contains vital information about the history of galactic evolution. One of the main 
problems associated with estimating an isotropic CIRB in the near infrared (1-5 //m) 
is the unknown contribution from stars within our own galaxy. The optimal obser- 
vational window to search for a background in the near-IR is at 3.5 //m since that is 
the wavelength region where the other main foreground, the zodiacal dust emission, is 
the least. It is not possible to map out the entire 3.5 /im sky at a resolution which 
will accurately estimate the flux from stars. However, since the CIRB is presumably 
isotropic, it can potentially be detected by selecting a smaller field and imaging it at 
good resolution to estimate the stellar intensity. We selected a 2°x2° "dark spot" 
near the North Galactic Pole which had the least intensity at 3.5 /xm after a zodiacal 
light model was subtracted from the all-sky maps generated by the Diffuse InfraRed 
Background Experiment (DIRBE). Still, the large area of the field made it very difficult 
to mosaic at 3.5 /xm using the available arrays. Thus, the field was mosaiced at 2.2 ^um, 
then the bright stars were selected and re-imaged at 2.2 and 3.5 (j.m.. The resulting total 
intensity of the bright stars was combined with a model for the contribution from dim- 
mer stars and subtracted from the zodi-subtracted DIRBE map. The contribution from 
the interstellar medium was also subtracted leaving a residual intensity at 2.2 //m of: 
16.4 ± 4.4 kJy sr"^ or 22.4 ± 6.0 nW m'^ sr-\ and at 3.5 /xm of: 12.8 ± 3.8 kJy sr'^ or 
11.0 lb 3.3 nW m~^ sr~^. The nature of our analysis suggests that this excess emission 
is probably a detection of the cosmic background in the near infrared. 

Subject headings: cosmology: observations — diffuse radiation — infrared: general 
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1. Introduction 

One of the few tools available for studying the "Dark Ages" following the era of matter- 
radiation decoupling is the Cosmic InfraRed Background (CIRB) radiation. The CIRB is the result 
of the cumulative, short wavelength emissions from pregalactic, protogalactic, and galactic systems 
which throTigh dust reprocessing and cosmological redshifting arc now at infrared wavelengths. 
Although the importance of looking for an extragalactic infrared background has been discussed 
for some time (Partridge and Peebles 1967; Low and Tucker 1968; Peebles 1969; Harwit 1970; 
Kaufman 1976), relatively little theoretical attention had been paid to it due to the difficulties 
inherent in trying to observe and verify theoretical predictions. But advances in infrared technology 
have stimulated an increasing amount of interest in trying to determine the characteristics of the 
CIRB (Bond et al. 1986; McDowell 1986; Fabbri et al. 1987; Fabbri 1988; Bond et al. 1991; 
Lonsdale 1995; Malkan &: Steckcr 1998; Pei et al. 1999). The theoretical work indicates that 
measuring the spectral intensity and the anisotropy of the CIRB will have important implications 
regarding the amount of matter undergoing luminous episodes in the pregalactic Universe, the 
nature and evolution of those luminosity sources, the nature and distribution of cosmic dust, and 
the density and luminosity evolution of infrared-bright galaxies (for review see Hauser (1996)). 

Although technology is no longer a hindrance in observing the CIRB, there are other difficulties 
in trying to see this cosmic relic. The bright foreground from the atmosphere of the earth, the dust 
in the solar system, and the stellar and interstellar emissions of our own galaxy combine to make 
detecting the CIRB a formidable task. 

Several experiments have been carried out to try to determine the CIRB from within the 

atmosphere. Matsumoto et al. (1988) and Noda et al. (1992) used rocket borne cameras to try to 
escape the contamination from the atmosphere and observe the CIRB. They had only limited success 
and could set upper limits on the CIRB without any detections. Lower limits were determined 
from deep galaxy counts conducted in the near-IR. These two limits began setting the first tight 
constraints on galaxy evolution models. 

The logical next step in trying to detect the CIRB was to go into space to try and eliminate the 
more local foreground; thus, the Cosmic Background Explorer (COBE) satellite carried with it an 
experiment specifically designed to detect the CIRB (Boggess et al. 1992). The Diffuse InfraRed 
Background Experiment (DIRBE) was cryogenically cooled, thus allowing it to observe multiple 
infrared wavelengths without much instrumental background. Having eliminated the eff"cct of the 
earth's atmosphere and the instrumental background, only the zodiacal dust in the solar system 
and the stellar and interstellar emissions of the Galaxy remained. 

Even with the contamination of the atmosphere removed, the DIRBE data did not readily 
reveal an isotropic background. Removing the remaining foregrounds turned out to be a difficult 
task with the zodiacal dust proving to be especially difficult as none of the models were able to 
completely eliminate it from the data, especially in the near and mid-infrared where the dust 
emits and scatters the most. The emissions from interstellar dust were more readily modeled 
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Fig. 1. — Limits and detections by DIRBE compared to various CIRB models. The DIRBE 2a 
upper limits at 1.25-60 ^m. are represented by solid downward triangles. The detections at 2.2 /xm 
and 3.5 jim. (this paper), and at 140 ^im. and 240 jim. (Hauser et al. 1998) are represented by solid 
circles with ±2c7 error bars. The 100 /xm intensity is represented by both the 95% confidence limit 
(open circle) on the non-galactic, non-zodiacal intensity of 5-35 nW m^^ sr^^ and the quoted upper 
limit on the cosmic background. The Far InfraRed Absolute Spectrophotometer (FIRAS) 125-5000 
/xm detection (Fixsen et al. 1998) is shown by a dotted curve while the Lagache et al. (1999) is 
shown by a dashed curve; the A < 2.2 fim lower limits of Pozzetti et al. (1998) are represented 
by small upward triangles. The upper limits on the CIRB based on TeV observations of Mrk 501 
(Stanev & Franceschini 1998) arc shown as small downward triangles, while the "detection" based 
on TeV observations of Mrk 421 (Dwek & Slavin 1994) is shown as a dashed line. The gray band 
shows the models of Malkan & Stecker (1998). 
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and removed and the emissions from stars were also modeled and removed. These attempts at 
foreground removal finally yielded values for a CIRB at 140 /im and 240 fim, first as a possible 
detection from an independent group (Schlegel et al. 1998) and then the definitive value from the 
DIRBE science team (Hauser et al. 1998). These values are shown in Figure 1 in relation to various 
models. As can be seen, the models are within factors of 2-6 of the observed values at 140 fim and 
240 /xm, indicating a convergence between observation and theory in the far-IR. 

Yet in the near-IR (NIR) (1-5 fim), DIRBE was unable to detect the CIRB. The contamination 
from the zodiacal dust as well as uncertainties in the contribution from galactic stars resulted in 
residuals which had large error bars and thus could only be considered as upper limits (Hauser et 
al. 1998). In addition, the residual maps in the NIR failed the tests for isotropy even in limited 
regions of the sky. 

1.1. Zodiacal Light Removal 

The zodiacal dust is made up of dust from asteroids and comets and extends out to the asteroid 
belt. The particles scatter and emit light from the UV through the infrared with a minimum of 
the sum of scattering and emission at 3.5 /xm. In recent decades a great deal has been added to 
our knowledge of the zodiacal dust, especially with the advent of satellites like IRAS and COBE. 
Taking the sum of our knowledge today, we know that the zodiacal dust is distributed in a fan 
shape centered on the sun and that it is slightly tilted with respect to the ecliptic. There are dust 
bands associated with astcroidal families, and there is a zodiacal dust ring at a distance of 1.01 AU. 
Unfortunately this accumulated knowledge has not brought with it a completely accurate model 
of the dust distribution. The DIRBE science team has produced a model of the dust distribution 
(Kelsall et al. 1998) to subtract from the DIRBE all-sky maps, but the model parameters are not 
unique and thus, after subtraction, the model zodiacal light leaves some residual effects. 

The Kelsall et al. (1998) model leaves a large residual intensity in the galactic polar caps at 
25 fim, the DIRBE band that is most dominated by the zodiacal light. For example, the 25 /xm 
intensity toward the DIRBE dark spot at {l,b) = (120.8°, 65.9°) in the DIRBE Zodi-Subtracted 
Mission Average (ZSMA) maps is 1.76 MJy sr^"*^. This cannot be a cosmic background because the 
lack of 7-ray emission toward Mkn 501 limits the CIRB to be < 33 kJy sr'^ (Funk et al. 1998). 
It also cannot be galactic cirrus because the 100 /tm intensity in this field is 1.27 MJy sr~^ in 
the ZSMA maps, and Arendt et al. (1998) specify the ISM intensity as i?(A)(7(100) - /o), with 
i?(25) = 0.0480 and = 0.66 MJy sr^^ so the ISM intensity is 29 kJy sr^^ By ehmination, most 
of this intensity must be zodiacal. 

In order to reduce the residual zodiacal emission in the maps, Wright (1997) added one "ob- 
servation" that the high b intensity at 25 /xm should be zero to the more than 10^ observations 
used in the zodiacal model fitting. Even this very low weight pseudo-observation lowered the 25 
/xm intensity in the dark spot to 0.26 MJy sr~^. This indicates that the isotropic component of the 
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zodiacal emission is very poorly constrained in fits that just look at the time variation to measure 
the zodiacal light. The Appendix in Wright (1998) discusses zodiacal light models in more depth, 
and the Appendix of this paper gives the actual parameters of the model we have used. 

Since the residual 25 ^um intensity is now only 1% of the total zodiacal emission, wc might 
hope for errors in the ZL model equal to 1% of the ecliptic pole intensity. But the situation is 
more uncertain at 2.2 and 3.5 jim due to the scattered component of the ZL. Adjusting the thermal 
emission component to fit the 25 fim intensity will not necessarily lead to a correct scattered 
component. So we have adopted ZL modeling errors of 5% of the intensity at the ecliptic poles at 
2.2 and 3.5 /xm. These errors are slightly lower at 2.2 //m (5.2 vs. 6 nW m~^ sr~^) and higher at 
3.5 /xm (2.8 vs. 2 nW m~^ sr~^) than the errors adopted by Kelsall et al. (1998). 

1.2. Interstellar Medium Contribution 

Another foreground is the emission from dust in the interstellar medium (ISM) . This contribu- 
tion is significant in the FIR but is relatively small in the NIR. However, it can not be completely 
neglected. The distribution of the dust was modeled by Arendt et al. (1998) by correlating 
the DIRBE 100 /xm map with external data sets of gas and dust like HI, HII, and CO. Arendt et 
al. (1998) specify that the ISM subtracted map at wavelength A is 1, 6)— i?(A)(7i,(100, 1, b)—Io). 
lo = 0.66 MJy sr~^ is the non-zodiacal, non-galactic component of the 100 /zm map. The coefficient 
R{2.2 iim) is zero, but i?(3.5 fim) = 0.00183 (Arendt et al. 1998) due to single photon induced 
thermal fluorescence in the 3.3 fim PAH feature in very small grains. We have made this correction 
on a pixel by pixel basis, and it should remove the effect of galactic cirrus if the global coefficient 
applies to our spot. The ISM correction listed in Table 2 is the correction averaged over the 17 
pixels in the "dark spot" . 

1.3. Foreground Starlight Removal 

Since there is no all-sky map at a sufficiently high resolution in the NIR, it is difficult to 
calculate the cumulative intensity of stars at varying galactic locations. This calculation is vital in 
eliminating one of the main foregrounds to the CIRB between 1 /xm and 5 /xm. 

To get an estimate of the intensity from stars, theoretical models of the Galaxy have been 
constructed from which number counts in any direction can be derived. The most recent and 
successful model of the Galaxy was an enhancement of the Wainscoat et al. (1992) model made 
by Cohen (1994) who divided the Galaxy into 5 main components: the disk, the spiral arms, the 
molecular ring, the central bulge, and the extended halo, as well as several minor components: 
Gould's belt, local molecular clouds, reflection nebulae, etc. In his model, each major component 
has 82 different stellar types distributed among the components with weighting techniques designed 
to reproduce observations of each component in the NIR. Additionally the spectral classes are 
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characterized by absolute magnitudes at IR wavelengths, a magnitude dispersion, a scale height 
above the galactic plane, and a solar neighborhood density. Finally, dust in the Galaxy is modeled 
using the Rieke & Lebofsky (1985) uniform extinction law. 

For the purpose of detecting the CIRB, the model has to be accurate at high galactic latitudes 
where the number counts are lowest and present the least amount of contamination. This is where 
the difficulties with the model become evident. Since the sun is situated about 15 pc above the disk 
and we are looking through the disk and the halo when we are looking at high galactic latitudes, 
any error in the halo-to-disk weighting will generate some error in the star counts. An additional 
source of error will be the actual distance of the sun above the disk. The greater the height above 
the disk, the lower the contribution of disk stars when looking North vs. South. These problems 
affect the counts and the surface brightness calculations at about the 10% level (Cohen 1996). 

A different source of error arises when the inherently statistical predictions of the model are 
subtracted from observed data. Any model will not exactly recreate the number counts in the 
galactic region of interest. There will be a variance from the actual number of stars. This variance 
is of little consequence for dim stars since each star contributes relatively little to the overall 
intensity, but the variance becomes very important for the brighter stars where the flux of each 
individual star contributes significantly to the overall intensity. 

While the zodiacal light model is intrinsically on the same flux scale as the DIRBE data, 
the relative calibration of the external stellar data for bright stars and models for faint stars 
vs. the DIRBE data needs to be known with good accuracy. The faint star model of Arendt et 
al. (1998) used fluxes at zero magnitude of Fo{K) = 612.3 Jy and Fo{L) = 285 Jy (COBE DIRBE 
Explanatory Supplement 1998). We have checked these values by extracting the total flux in 16 
pixel areas around the bright stars /3 And, a Tau, a Aur, a Ori, a Boo, a Her Sz (3 Peg. The 
median Fo{K) derived from these stars is Fo{K) = 614 Jy with a mean and standard deviation of 
the mean 616.6 it 12.5 Jy. This calibration agrees with Arendt et al. (1998) to better than 1%. But 
at 3.5 fiui, the median Fo{L) = 263 Jy and the mean is 264.9 ib 3.4 Jy. This is an 8% discrepancy. 
We will use our median value, Fo{L) = 263 Jy, and specify the changes created if the Arendt et 
al. (1998) value is used instead. 

Until recently, only this modeling approach was used to eliminate the galactic foreground. 
However, Dwek &; Arendt (1998; DA98) attempted to address the uncertainty in the number counts 
of stars using an innovative technique. They used the zodiacal light subtracted DIRBE all-sky maps 
at 2.2 /xm to help determine the background at 3.5 /xm. First they subtracted a background from 
the 2.2 /xm data based on deep K band (2.2 /xm) galaxy counts. This then should have left only 
the stellar contribution in the 2.2 /xm maps. Then the 1.25, 3.5, and 4.9 /xm maps were used to 
determine the proper colors of the stars which were then subtracted from the 3.5 ;um map. Based 
on this approach they made a tentative claim of detection for the CIRB at 3.5 /xm of: 

1^/^(3.5 /tm) = 9.9 + 0.312(z//^(2.2 /xm) - 7.4) ± 2.9 (1) 

where uli, is in units of nW m~^ sr~^. The factor of 0.312 is the color correction which was found 
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from the correlation of the 2.2 fim map with the 3.5 /um map, and the (1^/^(2.2 /xm) — 7.4) is the 
real value of the CIRB at K minus the assumed value of the CIRB of 7.4 nW m~^ sr~^ from the 
deep galaxy counts. 

1.4. Motivation for this Project 

As mentioned earlier, the strongest constraints available on the CIRB in the NIR are the upper 
limits by the DIRBE science team (Hauser et al. 1998), and the tentative detection at 3.5 //m by 
DA98 also based on DIRBE data. The DIRBE team did not make a detection claim because 
they did not have a significant residual signal above the uncertainties associated with the stellar 
subtraction and also because they had not achieved isotropy in their residual maps (Arendt et 
al. 1998). Their stellar subtraction was achieved in 2 steps. In the first step they removed all stars 
brighter than 15 Jy since they were bright enough to be identifiable in the low resolution DIRBE 
image. This translated into a brightness cutoff at K of 4*^ magnitude and at L of 3^'^ magnitude. In 
the second step they modeled the contribution of stars dimmer than the cutoff value using Cohen's 
model of the Galaxy. Our intention was to reduce those modeling uncertainties by imaging one 
section of sky at high resolution to a limit at least 2 orders of magnitude fainter than the Arendt 
et al. (1998) bright star cutoff. This translated into 9*^ magnitude at K and 8*'' magnitude at L. 
At these limits, the standard deviation in the intensity contributed by the Poisson fluctuations in 
the numbers of faint stars is reduced to a level 16 times smaller than the standard deviation for 
the Arendt et al. (1998) threshold. 

DA98 claimed a tentative detection at 3.5 //m. Although they seemed to have isotropy, they 
did have a low-level large scale gradient, probably due to zodiacal light, that remained in their 
images. In addition, their value at 3.5 /xm was still based on an assumed value of the CIRB at 2.2 
/xm obtained from ii'— band galaxy counts. 

Our intention here was to use a much more stringent zodiacal light subtraction routine as well 
as to avoid any assumption about the CIRB at other wavelengths. With these two new approaches 
in mind we hoped to place a much tighter constraint on the CIRB at 2.2 and 3.5 /xm. 

2. Observations 

2.1. The DIRBE Data 

The primary dataset for this project was the all-sky map produced by the DIRBE instrument 
on the COBE satellite. DIRBE was helium cooled and observed the sky at 10 different IR wavebands 
from 1.25 /im to 240 /vm. The sky was mapped onto 393,216 pixels with a pixel scale of 0.32° pixel"^ 
and each pixel was observed several hundred times. The instrument was also equipped with a 
chopper which chopped between the sky and an internal calibrator 32 times per second. Finally, the 
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Fig. 2. — All-sky map from DIRBE at 3.5 fj,m. in galactic coordinates with intensities in MJy sr~^. 
The S-shape of the zodiacal light is evident as well as the contamination by individual stars which 
appear almost as noise on the images. The "dark spot" (a = 12^ 56.9™, 5 = 51° 11' (J2000)) and 
the region used for the 2MASS analysis (a = 13^ 24™, (5 = 15° (J2000)) are outlined. 

instrument was re-calibrated 5 times per orbit. All this redundancy resulted in the best absolutely 
calibrated map of the infrared sky obtained to date (Figure 2). 

The first step in determining the CIRB from the DIRBE data was to determine the location 
in the sky with the lowest intensity. But since the Earth sits within the zodiacal dust cloud, the 
region of lowest emission was still heavily contaminated by scattering and emission from that dust. 
Figure 2 shows the prominent S-shape (in galactic coordinates) of the zodiacal dust, and before 
anything else could be attempted with the data, that contamination needed to be removed. To 
remove the contribution of the zodiacal light required modeling the distribution of the dust as 
well as its emission and scattering characteristics. This task was accomplished by Wright (1998) 
resulting in maps where the remaining contamination was from the Galaxy. We did not use the 
zodi-subtr acted maps generated by the DIRBE science team (Kelsall et al. 1998) because this 
model leaves a large residual 25 /xm intensity in the region of our ground-based survey. Thus we 
used a model constructed using the "very strong no-zodi" condition defined by Wright (1997), 
which reduces the residual 25 /xm intensity by a factor of 7. 

From the cleaned map, the area of lowest emission, a "dark spot" , was chosen which was a 
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region spanning 2 degrees by 2 degrees. This 2°x2° area of lowest emission turned out to be near 
the North Galactic Pole at a = 12^ 56.9"", and 5 = 51° 11' (J2000) {I = 120.8°, b = 65.9°). The 
low value at this location represents a combination of the CIRB and the unresolved stars within the 
large (0.7°) DIRBE beam. The next step was then to determine the contribution from unresolved 
stars using higher resolution ground based imaging. 

2.2. The Initial Map 

Obtaining the fluxes of stars in the determined "dark spot" at 3.5 /xm was the primary goal of 
this project, but mosaicing an area of 4 square degrees is not possible with present day technology. 
Modern arrays do not read out fast enough to be able to image a large region of sky before the 
background at 3.5 /7,m saturates the array; therefore, most cameras that are meant to image in the 
thermal infrared (> 3 fim) are designed to either image a very small region of sky to minimize the 
background or to read-out only a small subarray of the detector chip. 

The main camera available for this project was the UCLA Twin Channel Infrared Camera 
(known locally as Gemini) (McLean et al. 1993) which mounts onto the Lick 3-m telescope. 
On the 3-m, Gemini's two 256x256 arrays have identical 3'x3' fields of view, one array covering 
the wavelength range from 1 to 2.5 /xm with a HgCdTe array, and the other array covering the 
wavelength range from 1 to 5 /xm with an InSb array. The camera employs the subarray technique 
for readout when imaging at wavelengths longer than 3 /im, outputting 64x64 pixels giving it a 
field of view of only 45"x45". With a field of view this small it would be impossible to mosaic the 
required area of the "dark spot" which covers the equivalent area of 20 full moons. 

The solution that was employed was to image the area at K with a different camera which due 
to the lower sky background at K, could have a much larger field of view, then choose the brightest 
stars from that sample to re-image at 3.5 /xm with Gemini on the 3-m. The camera available for this 
approach was the Lick InfraRed Camera 2 (LIRC2). LIRC2 could mount onto the 1-m telescope 
at Lick and have a field of view of 7.29' x 7.29', easily allowing the 4 square degrees to be mapped 
in a short period of time. This was done over the nights of 1997, March 28 and 29. The 302 stars 
catalogued in the LIRC2 survey are shown in Figure 3. 

The process was complicated during the reduction of the LIRC2 data when it was discovered 
that LIRC2 had photometric problems which precluded highly accurate photometry. The accuracy 
was still sufficient to estimate relative brightnesses, which was the main result that was required of 
the data. Having obtained the list of stars from the brightest to the dimmest in the "dark spot" 
we could now proceed with imaging them in the L band. 
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2.3. The Main NIR data 

We had three clear and photometric nights at the Lick 3-m with Gemini on the nights of 1999, 
March 26 to 28. Using Gemini's twin arrays we imaged our Ust of stars in both the K band and 
the L band simultaneously. 

The data were reduced using standard data reduction routines with the IRAF^ software pack- 
age. Since the sky brightness varies dramatically over very short periods of time at 3.5 /xm, differ- 
encing of sequential images was used in the data reduction to remove detector signature and sky 
background. Each star was imaged at 4 different locations on the array, while the 10 brightest stars 
were imaged a total of 12 times. 

We also obtained time at the 1.5-m telescope at Palomar on the night of 1999, April 4, using 
its IR camera. Wc imaged the 20 brightest stars at K band as a check of our Gemini data as well 
as imaging the list at J band. These data were reduced in a similar fashion to the Gemini data. 
The uncertainties for both data sets were determined by the standard deviation of the brightness 
of the standard star throughout the night. 

A total of 77 stars out of the original sample of 302 were observed either at Palomar or with 
Gemini or both. We have restricted our analysis to the brightest part of the sample in order to 
have a reliable catalog with a well defined limit. Table 1 shows the 28 bright stars with m < 9 in 
the dark spot, with the original LIRC2 K magnitude, the Palomar J and K magnitudes (J60 and 
K60) and the Lick Gemini K and L magnitudes. We have used the Lick K and L magnitudes in 
this paper. 

3. Combining the two Datasets 

Having obtained the photometry from the ground for the bright stars, the task was to match 
up the ground based data to the space based data. Figure 3 shows a projection of the DIRBE 
pixels on the "dark spot" as well as the stars imaged by LIRC2. 

The most obvious difficulty turns out to be that the DIRBE pixels are not aligned with the EW 
and NS axes; therefore, some only partially enter the 2°x2° box. Also the values within the pixels 
themselves represent the average intensity seen when the 0.7° square DIRBE beam is centered 
within the pixel boundary, and hence are affected by the pixels around them since the beam does 
not abruptly stop at the edge of the pixel. 

To compute the response in a given pixel to a given star, we found the probability that the star 



^IRAF ( Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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Fig. 3. — A superposition of the stars within the surveyed 2°x2° area (thin-hned box gray box) 
and the DIRBE pixels (thick-hncd skew rectangles). The 17 DIRBE pixels used in our analysis are 
shaded in darker gray. The numbers in the pixels are from the top: The DIRBE pixel number and 
the J, K, and L intensities after zodi subtraction in kJy sr~^. The grid lines are at even degrees of 
a and 5. The circles are the stars imaged with LIRC2. The bigger the size, the higher the flux. 
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was inside the DIRBE 0.7° square beam assuming that the beam center was uniformly distributed 
inside the pixel, and that the beam orientation was uniform in angle. This is a fair assumption 
since each pixel was observed hundreds of times in many different angles during the lifetime of the 
mission. To be specific, let H(b, 6, s) be 1 if the star with position given by the unit vector s is in 
the square beam when it is centered at position b and oriented with position angle 9, and if the 
star is not in the beam. Then the probability of the j*^ star being in data taken in the z*'* pixel is 
given by 

where is the angular extent of the i^^ pixel. Then the predicted contribution from bright stars 
to the intensity in the i*'* pixel is given by 

Bi = n-^Y.P^,F, (3) 
j 

where Jl is the solid angle of the DIRBE beam and Fj is the flux of the j*'* star. We also need 
to know the variance of the intensity in the pixel due to the random orientations and random 
placements of observations within the pixel boundaries. This can be obtained from the expected 
value of within a pixel, but since = H we find that the variance of pij is given by Pij{^—Pij)- 
Thus stars at the edge of the beam (p ~ 0.5) make the greatest contribution to the variance, a 
phenomenon known as "flicker noise" . We then get 

a{Bif = n-^J2^ij{l-pij) +p%e^]Fj (4) 

3 

where e is an uncertainty in cross-calibrating the DIRBE and ground-based data for individual 
stars which also includes an allowance for variability in the stars. We use = 0.1 in this work. 
Many IR bright stars are long period variables and will not have the same flux during the DIRBE 
observations as they have nine years later during our ground-based survey, so a fairly large value of 
is appropriate. Note that the variance sum is dominated by the brightest stars, and is minimized 
in the "dark spot" . 

We have evaluated the bright star contribution for all stars with < 9 mag or L < 9 mag 
in the dark spot. Note that the cutoff magnitude does not come from LIRC2 so as to avoid any 
additional error being introduced because of LIRC2's photometric inaccuracy in the final result. 
Stars were detected below these limits with accurate photometry but were not used, with the 
conservative cut-off being at 9th magnitude. There are only 17 DIRBE pixels whose response to 
a uniform illumination pattern is at least 90% due to radiation from within the dark spot, and 
these 17 pixels form our DIRBE sample. Given the DIRBE data Di and the bright star intensity 
5j, we can estimate the cosmic infrared background using Ci = Di — Bi — F — Z , where F is the 
contribution from faint stars evaluated from star count models, and Z is the zodiacal contribution. 
Both F and Z are essentially constant over our spot and thus are not subscripted by pixel. We 
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Fig. 4. — The zodi-subtracted DIRBE 3.5 /xm intensity plotted vs. the predicted intensity from 
bright stars based on Lick data for stars with L < 9 mag in the dark spot at b) = (120.8°, 65.9°). 
The hne is of unit slope. 
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Fig. 5. — The zodi-subtracted DIRBE 2.2 fiui intensity plotted vs. the predicted intensity from 
bright stars based on Lick data for stars with K < 9 mag in the dark spot at b) = (120.8°, 65.9°). 
The hne is of unit slope. 
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2MASS intensity [l^Jy/sr] 

Fig. 6. — The zodi-subtracted DIRBE 2.2 /xm intensity plotted vs. the predicted intensity from 
bright stars based on 2MASS stars with K <9 mag in the region at (l,b) = (337.1°, 75.6°). The 
arrow labeled K=5 shows the intensity produced by a single 5*'* magnitude star in the DIRBE 
beam. The cross-hatched area in the lower left shows the plotting box for Figure 5. The line is of 
unit slope and the intercept is in good agreement with our value for the intensity from stars K > 9 
mag plus the CIRB. 
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combine the several estimates of the cosmic infrared background using a weighted median: we find 
the value C that minimizes 

Cj-C 
^-^^ 

The uncertainty in C is found using the half-sample method (Babu &: Fcigclson 1996): many Monte 
Carlo runs were performed each with a randomly chosen half of the pixels, and the scatter in these 
determinations gives the uncertainty in C. The result for the CIRB plus the faint stars {C+F) at K 
is 28.27ib0.90 kJy sr~^, and at L is 17.81ib0.75 kJy sr~^. These uncertainties are based only on the 
statistical pixel to pixel scatter within the "dark spot" , but they do include the fluctuations in the 
faint star contribution. Figure 4 shows the zodi-subtracted DIRBE 3.5 //m intensity, Di — Z, plotted 
against the bright star intensity predicted from the Lick L band data. The slope of a line fitted to 
these points is only 0.4, but this does not mean that the DIRBE calibration differs by a factor of 0.4 
from the ground-based calibration. Rather it shows that the variance of the stellar contribution is 
not dominated by the bright stars in this "dark spot" , and thus the predicted stellar intensity is a 
noisy estimator of the actual stellar intensity. In this circumstance, the slope of a least-squares line 
is biased low, as shown in Figure 2 of Wright (1992). Figure 5 shows the zodi-subtracted DIRBE 
2.2 nm intensity vs. the bright star intensity predicted from the Lick K band data. Figure 6 shows 
a similar plot using 2MASS data in a 2.4° x 4° region centered at a = 13'' 24"'', 6 = 15°. This region 
is at high galactic latitude (Z = 337°, b = 76°) but has not been selected to have very few bright 
stars. Even though pixels affected by stars listed as saturated in the 2MASS catalog have been 
deleted, there appear to l)c several stars that are partially saturated in the 2MASS data leading 
to pixels with high DIRBE intensities but low 2MASS intensities. However, the plotted line with 
unit slope is clearly consistent with a large number of pixels. The intercept is 29.9 ± 5.9 kJy sr~^ 
which compares very well with the 28.3 it 0.9 kJy sr~^ found in the "dark spot". 

The 17 pixels are not independent because a given star affects several pixels. The ratio of the 
0.7° square beam solid angle to the pixel solid angle is 4.7:1. But the effective beam has fuzzy 
edges where Pij is close to 0.5, and the appropriate beam solid angle to use is the noise effective 
solid angle Jig = (YliPij)'^ / J2iPij = ^ pixels (Wright 1985). We therefore multiply the statistical 
uncertainty given above by a factor of a/6 ~ 2.5, giving the error bars on the m < 9 contribution 
listed in Table 2. But this uncertainty also includes the s/N fluctuations in the faint stars and any 
DIRBE instrument noise. 

We have run the Wainscoat et al. (1992) star count model for the "dark spot" and our 2MASS 
region. In the 2MASS region the model predicts 785 stars with 9 < K <\2 mag and we find 704 
stars. Thus the model count is 11.5% ± 3.8% too high. The predicted integrated intensity from 
stars with 9 < < 12 mag is 11.2% ± 5.4% too high. This 11% difference is well within the 10-15% 
estimated uncertainty in the model (Arendt et al. 1998). In the absence of data similar to 2MASS 
at 3.5 /Ltm, we assume that the model colors are correct, and apply this correction factor to the 
model values for the contribution from stars fainter than 9*^ magnitude in the "dark spot" at both 
2.2 /xm and 3.5 /xm and get F = 11.9 ± 0.6 kJy si'^ at K and F = 5.7 ± 0.3 kJy sr'^ at L. 



(5) 
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There is a bias in F introduced by our selection of a "dark spot" . Let the total flux which we 

select on he T = B + F , where the bright star contribution B and the faint star contribution F are 
independent random variables. If B and F are Gaussian, and we select on T, then there is a bias 
in F oi ap X (o'f/c^t) x (T/ctt)- We have selected a dark spot that is a fluctuation with probability 
of ~ 10"^. The actual fluctuations are not Gaussian, and B has infinite variance, but if we limit T 
to stars fainter than the 15 Jy cutoff used by Arendt et al. (1998) a Gaussian approximation is not 
too bad. Then P = 10~^ is a — 4cj selection. At 2.2 fim, ax = 10.8 kJy sr~^ and ap = 0.66 kJy sr"-*- 
in a 2° X 2° box, leading to a bias in F of 0.66 x (0.66/10.8) x (-4) = -0.16 kJy sr'^ At 3.5 ^m, 
Gt = 8.8 kJy sr~^ and ap = 0.32 kJy sr~^, leading to a bias in F of —0.05 kJy sr~^. We have not 
corrected for these biases, but in principle our CIRB results should be increased by 0.16 kJy sr~-^ 
at 2.2 jjim and 0.05 kJy sr~^ at 3.5 /xm. 

There is an ISM component at 3.5 jim. that also needs to be subtracted. It is 0.00183 times the 
difference between the zodi-subtracted 100 /itm intensity, which has a mean of Jioo = 1-02 MJy sr"-*^ 
in our "dark spot", and the zero Nh intercept of the 7ioo vs. Nh correlation, which is lo = 
0.66 MJy sr~^ (Arendt et al. 1998). Our maps used for the pixel by pixel bright star fitting 
had 0.00183 times the entire zodi-subtracted 100 /im intensity subtracted without first removing 
Iq. Furthermore, Iq using our zodi model would be 0.4 MJy sr~^ since we subtracted more zodi. 
Since we used 13.5 — 0.00183/ioo but want 13.5 — 0.00183(/ioo — -^o), we need to add a constant 
O.OOI83/0 = 0.7 kJy sr~^ to our 3.5 /xm results to complete the ISM correction. 

Noise due to the DIRBE detectors would be included in the pixel to pixel scatter, but inde- 
pendent from pixel to pixel. Thus the contribution from DIRBE noise should not be multiplied 
by the factor of 2.5 discussed above, and we have been conservative by assuming that the DIRBE 
noise is zero. There is no calibration error between DIRBE and the zodiacal light models since the 
models are fit to data on the DIRBE scale. 

Other errors could potentially arise from the zodiacal modeling uncertainties and the cross 
calibration uncertainties between the ground based and space based data. These uncertainties are 
dominated by errors in the zodiacal cloud model, and we estimate them to be 3.8 kJy sr^^ at 
2.2 /Ltm and 3.3 kJy sr~^ at 3.5 fim. Taking the quadrature sum of 2.5 times the statistical error 
and our estimated zodiacal modeling errors gives a net uncertainty of 4.4 kJy sr~^ at 2.2 /xm and 
3.8 kJy sr~^ at 3.5 /xm. 

Having determined the stellar contribution we removed it from the DIRBE map arriving at a 
residual signal at K of: 

I2.2^.m = 16.4 ± 4.4 kJysr-1 (6) 

or 

i^h.2tim = 22.4 ± 6.0 nW m-^ sr-^ (7) 

and a signal at L of: 

h.5^lm = 12.8 ± 3.8 kJy sr-i (8) 
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or 

r^h.bum = 11.0 ± 3.3 nW m-^ sv'^ (9) 

Table 2 shows the total DIRBE intensities and the zodiacal light, interstellar medium, bright stars, 
faint stars, and the residual signal. It should be emphasized that using the Kelsall et al. (1998) 
zodiacal light model instead of our model gives a residual intensity that is 5.1 kJy sr~^ higher at 
3.5 /xm and 6.6 kJy sr^^ higher at 2.2 ^um. 

Using the DIRBE faint source model calibration of Fo{L) = 285 Jy instead of the Fo{L) = 
263 Jy we measured would reduce the derived residual intensity at 3.5 jim by 0.74 kJy sr"-*-. 



4. Discussion 

The technique of zodiacal light subtraction used in this paper reduces the possibility that 
this residual signal is an artifact of this model. Certainly replacing our zodi model by the Kelsall 
et al. (1998) model would increase the apparent statistical significance of our result, but would 
leave an unexplained high latitude residual in the 25 /xm map. Our mapping of the dark spot also 

eliminates the possibility that the signal is stellar in origin: by observing and directly subtracting 
stars down to 9*^ magnitude we reduce our dependence on the Wainscoat et al. (1992) starcount 
model by a factor of 4 relative to the analysis of Arendt et al. (1998) and Hauser et al. (1998). 
Hence we conclude that the residual signal is the CIRB. 

Since we have only observed 0.01% of the sky, we cannot directly verify the isotropy of the 
CIRB that we see. But the evidence we do have for isotropy at 3.5 jum is qualitatively the same 
kind of evidence that Hauser et al. (1998) had for isotropy of CIRB at 240 /xm. The 240 /xm map 
is iVOT isotropic, but /240 — -R(240)(/ioo — -^o) is. Finding the level of this isotropic signal requires a 
determination of Jo by carefully subtracting the galactic signal in a "dark spot" : the Lockman hole. 
DA98 show that 13.5 — 0.312(/2.2 — -^o) is isotropic except for a low-level large-scale gradient caused 
by an error in the zodiacal light model. DA98 did not give the magnitude of the zodiacal gradient 
that they saw at 3.5 /xm, but Wright & Reese (2000) find the slope vs. esc (3 to be 3.1 nW m~^ sr"-*^ 
in the ZSMA maps used by DA98 and 2.3 nW m~^ sr~^ in the maps used in this paper. This 3.5 
/Ltm slope is less than the 3.6 nW m^^ sr"^ slope in the 240 /xm ISM2 map reported on line 15 of 
Table 2 in Hauser et al. (1998), and much less than the 6.9 nW m~^ sr~^ standard deviation of 
the 240 fiui ISM2 slope. Thus the bolometer noise in the 240 /xm ISM2 map could mask a larger 
gradient than the gradient seen in the low noise 3.5 //m map. The slope in the 3.5 //m map is 
21% of the CIRB, while the reported 240 /xm ISM2 slope is 26 ± 51 % of the CIRB. Furthermore, 
the 2 point correlation function of the 3.5 /itm residual (Dwek Sz Arendt 1999), as a fraction of 
the CIRB, is smaller than the uncertainty in the 2 point correlation function of the 240 /im. ISM2 
map. While DA98 obtained a lower limit on /q using galaxy counts at 2.2 /im, we have determined 
the value of lo by carefully subtracting the galactic signal in our dark spot. Thus the Hauser et 
al. (1998) determination of the far-infrared CIRB and our determination of the near-infrared CIRB 
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are logically equivalent. But one has more confidence in the Hauser et al. (1998) result because 

Arendt et al. (1998) were able to determine their Iq at 100 /um in 3 areas with consistent results, 
and because the zodiacal light (and its modeling error) are a much smaller fraction of the CIRB at 
240 iiva. than at 3.5 /um. 

Although the uncertainty for the 2MASS spot is large, we find that the value of the residual 
signal derived from the "dark spot" agrees with the 2MASS spot. In the future, we will be able to 
further verify the isotropic nature of this signal. Since 2MASS saturates on stars that are quite faint 
for DIRBE, 2MASS data can only be used in regions with no bright stars. The current 2MASS data 
release does not cover any of our DIRBE dark spots. However, future data releases will eventually 
cover the whole sky, allowing a good isotropy check on the CIRB at 2.2 /xm. 

Our accepted proposal to survey the darkest spot in the Southern sky at the South Pole 
SPIREX/ABU facility in the 3.5 //m band during the 1999 austral winter was not executed due to 
equipment problems. Surveying other dark spots at 3.5 /xm would be very valuable but also very 
difficult with ground-based telescopes. 

Even though our evidence for the isotropy of the near infrared background is not strong, the 
uncertainty in our final result is still dominated by the uncertainties in modeling the zodiacal light. 
Our CIRB values at both 2.2 and 3.5 /xm are 19% of the zodiacal light at the ecliptic poles, while 
the 240 /xm background found by Hauser et al. (1998) is 3.6 times larger than the zodiacal light 
at the ecliptic poles. Thus the near IR background is much more sensitive to zodiacal light errors 
than the far IR background. 

5. Conclusion 

We compare our values to the values obtained by Hauser et al. (1998). and DA98 in Figure 
7. Hauser et al. (1998) made no claim for the CIRB since they did not get isotropy over a 
large part of the sky. However they did publish values for regions of the sky where all foreground 
contaminations were expected to be lowest. They labeled them as high quality (HQ) regions. HQA 
covered approximately 20% of the sky and HQB covered approximately 2% of the sky. Their 
HQA value at K was 9.5 ± 8.8 kJy sr^^ (13 ± 12 nW m-^ sr-^) while their HQB values was 
10.9 lb 8.8 kJy sr~^ (14.9 ± 12 nW m^^ sr^^), where the uncertainties are the systematic errors. 
Our value falls within one sigma of both their HQA and HQB values though that is due mostly to 
the fact that their error bars are very large. 

At L band their HQA value was 13±7kJy sr-i (ll±6nW m ^ sr ^) and their HQB value 
was 13.3 ± 7 kJy sr-i (11.4 ± 6 nW m ^ sr ^). Our L band value is also within one sigma though 
much closer to their quoted numbers. 

Since DA98 were able to make a much more precise measure than Hauser et al. (1998), the 
comparison between our value and theirs is very interesting. Inserting our value for i'Iy(2.2 nm) 
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Fig. 7. — Comparison of our results (point with error bars) and previous limits: the hatched region 
shows the Hauser et al. (1998) upper limits, while the gray bar shows the DA98 correlation. The 
point in the upper right shows the actual average DIRBE intensity in our 17 pixel dark spot, with 
the dashed line showing the zodiacal light subtraction, and the solid line showing the star light 
subtraction. 
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into the DA98 correlation predicts ul,y{3.5 /xm) = (14.6 it 3.5) nW sr~^ which is one standard 
deviation higher than our result. Note that DA98 essentially only determined that the high galactic 
latitude intensity could not be entirely due to galactic stars because a smooth component with 
a redder K-L color was needed. In principle, this "horse of a different color" could have been 
incompletely subtracted zodiacal light, since the K-L color of the zodiacal light is much redder 
than typical stars due to the contribution of thermal emission to the 3.5 /xm intensity. If we had 
used the Kelsall et al. (1998) zodiacal light model instead of our own, the inferred 2.2 //m CIRB 
would have increased by 6.6 kJy sr^^and the 3.5 /um CIRB would have increased by 5.1 kJy sr~^. 
These changes would have induced a change of 3 x (5.1/3.5 — 0.312 x 6.6/2.2) = 1.6 nW m^^ sr~-^ 
in the difference between the DA98 result and our result, in the sense that the difference would 
go from (3.6 it 3.5) nW m~^ sr~^ to (2.0 it 3.5) nW m~^ sr~^. Thus our values are in very good 
agreement with the DA98 correlation. Since we use a different zodiacal light model and a different 
method of star subtraction, the fact that we and DA98 have arrived at similar numbers gives 
confidence that the residual signal we derive at 2.2 ^um and 3.5 /Ltm is very likely the CIRB rather 
than an unknown component of Galactic emission. 

The COBE datascts were developed by the NASA Goddard Space Flight Center under the 
guidance of the COBE Science Working Group and were provided by the NSSDC. This publication 
makes use of data products from the Two Micron All Sky Survey, which is a joint project of 
the University of Massachusetts and the Infrared Processing and Analysis Center, funded by the 
National Aeronautics and Space Administration and the National Science Foundation. 

A. Zodiacal Light Model 

The zodiacal light model used in this paper is derived from the model described in Wright 
(1998). Three changes have been made: 

1. The scattering phase function coefficients p2o and p2i are allowed to be different for each of 
the three short wavelength bands. In Table 3, I use JKL to distinguish the values for the 
1.25, 2.2 and 3.5 /xm bands. 

2. The scaling of the IRAS zodiacal bands, pi^, is allowed to be different for the scattered light, 
denoted P135, than it is for the thermal emission, denoted pisr- 

3. The model was fit to the "very strong no zodi" condition of Wright (1997). 

The emissivities are given in Table 4 and the IRAS band parameters are given in Table 5. 
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Table 1. Photometry of Stars in the Dark Spot with K <9 mag or L < 9 mag 



Name 


Q!2000 


<^2000 


LIRC2 


J60 


K60 


K 


L 


D48 


13^02™53.7'* 


51°02'21" 


7.31 


6.69 


6.02 


6.06 


5.83 


C61 


12^56"'31.5'' 


51°04'49" 


6.36 


7.11 


6.36 


6.42 


6.17 


A96 


13'*03"*09.P 


51°42'28" 


6.83 


7.19 


6.55 


6.65 


6.48 


C250 


12'*52"*21.4^ 


50°16'56" 


7.63 


7.78 


7.51 


7.58 


7.51 


D240 


12'»59"*31.9^ 


50°21'47" 


8.39 


8.89 


8.22 


8.22 


7.67 


B192 


12'^54™45.0^ 


51°33'50" 


8.08 


8.43 


7.89 


7.97 


7.72 


A207 


12'^ 58"" 17.5'' 


51°16'28" 


8.11 


8.82 


8.04 


8.13 


7.72 


B262 


12'*52'"47.8^ 


51°12'09" 


8.14 


8.56 


8.29 


8.32 


7.84 


A210 


12'*58"*17.5^ 


51°16'28" 


9.25 


8.84 




7.42 


7.87 


D78 


12'*57""30.5^ 


50°49'26" 


8.10 


8.47 


7.96 


8.04 


7.88 


D74 


12'^56™52.3^ 


50°48'05" 


7.72 


8.45 


7.86 


7.94 


7.89 


A13 


13^00"00.2*' 


52°08'03" 


8.45 


9.49 


8.82 


8.85 


8.11 


C199 


12^50"56.4*' 


50°29'56" 


8.02 


8.75 


8.41 


8.48 


8.32 


A57 


13^01""21.r 


51°51'50" 


8.82 


9.55 


8.79 


8.87 


8.38 


A107 


12^59^440^ 


51°42'11" 


8.84 


9.94 




8.66 


8.42 


B257 


12"52"*25.8* 


51°14'13" 


9.02 


9.41 




9.04 


8.44 


B254 


12'*51""50.4'' 


51°16'29" 


8.59 


8.83 


8.52 


8.52 


8.46 


A86 


13^02'"56.6'' 


51°48'14" 


8.83 


9.89 




8.59 


8.47 


B267 


12^54'"52.6'' 


51°13'04" 


9.16 


9.92 




9.74 


8.50 


C261 


12^54™24.7*' 


50°15'43" 


8.80 


9.32 


8.75 


8.81 


8.72 


B50 


12'*56""43.5^ 


52°03'08" 


8.62 


8.92 


8.62 


8.62 


8.78 


A175 


13''00""40.4^ 


51°24'27" 


9.00 


9.23 




8.97 


8.88 


A117 


12'*57""52.r 


51°37'56" 


9.14 


9.71 




9.21 


8.90 


C211 


12^53™28.5'* 


50°28'24" 


8.86 


10.39 




8.82 


8.91 


C49 


12^52™44.3'' 


51°03'23" 


8.87 


10.37 




9.02 


8.93 


D221 


13''02'"37.7'' 


50°24'44" 


8.62 


9.31 


8.88 


8.88 


9.04 


B91 


12'*51""47.5^ 


51°51'15" 


8.28 


8.44 


8.16 


8.19 


9.07 


D60 


13'*01""03.8^ 


50°59'28" 


8.90 


9.13 




8.79 


9.51 
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Table 2. Decomposition of the DIRBE Intensity 



Component 


2.2 /xm 


3.5 /[/m 




kJy sr~^ 


kJy sr~^ 


Total 


137.5 ± 0.3 


105.3 ±0.3 


Zodi 


101.8 ±3.8 


80.4 ±3.3 


ISM 




1.1 ±0.2 


Stars, m < 9 mag 


7.4 ±2.2 


5.3 ± 1.8 


Stars, m > 9 mag 


11.9 ±0.6 


5.7 ±0.3 


CIRB 


16.4 ±4.4 


12.8 ±3.8 
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Table 3. Diffuse Cloud Parameters 



Parameter 


Value 


Description 


Pi 


1.2346 


radial density exponent 


P2 


0.4246 


radial temperature exponent 


PS 


3.5785 


vertical "scale height" 


P4 


0.9450 


vertical density exponent 


P5 


-1.3559 


In(sini) at break 


P6 


0.3838 


10 X cloud pole X component 


P7 


-0.0758 


10 X cloud pole y component 


P8 


-0.0195 


10 X cloud offset x component 


P9 


-0.0471 


10 X cloud offset y component 


PlO 


0.6098 


10 X density contrast of Dcrmott ring 


Pll 


5.5405 


In(TB), band temperature at R = 1 


Pus 


0.5722 


scattering scale factor for bands 


Put 


l.oooi 


thermal scale factor for bands 


Pl3 


0.3161 


the "dimple" in Dermott's ring 


pu 


7.8852 


vertical scale for Dermott's ring 


Pl5 


-0.0226 


spherical term in vertical density 


Pie 


0.0289 


sin^ i term in vertical density 


Pl7 


-0.0262 


additional density at sini 0.5 


Pl8 


-0.1977 


additional density at | sini| « 0.25 


Pl9 


-0.0294 


additional density at | sini| 0.17 


P2QJ 


-0.2616 


phase function linear coefficient at 1.25 /um 


P20K 


-0.2865 


phase function linear coefficient at 2.2 /Ltm 


P20L 


-0.2634 


phase function linear coefficient at 3.5 /im 


P21J 


0.6761 


phase function quadratic coefficient at 1.25 //m 


P21K 


0.5318 


phase function quadratic coefficient at 2.2 /xm 


P21L 


0.3988 


phase function quadratic coefficient at 3.5 fim 



Table 4. Scattering and Emission Efficiencies 



Parameter Value Description 



P22 


0.6647 


0"! scaling 


P23 


0.7762 


(72 scaling 


P24 


0.7499 


(73 scaling 


P25 


1.8247 


K2 scaling 


P26 


1.3150 


K4 scaling 


P27 


1.1088 


K5 scaling 


P28 


1.2525 


kq scaling 


P29 


0.8657 


Kr scaling 


P30 


0.7547 


Kg scaling 


P31 


0.8147 


kq scaling 


P32 


0.6130 


Kio scaling 



Tabic 5. IRAS Band Parameters 



Parameter 


Value 


Description 


qi 


1.7058 


10 X {smi)max for Band 1 


92 


0.1963 


Band 1 normalization 


93 


0.2722 


10 X {smi)max for Band 2 


94 


0.2515 


Band 2 normalization 


95 


0.1167 


10 X band pole x component 


96 


-0.2691 


10 X band pole y component 


97 


-0.9321 


10 X band offset x component 


98 


1.9164 


10 X band offset y component 


Ri 


3.14 


Outer radius for Band 1 


R2 


3.02 


Outer radius for Band 2 



